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Abstract 

We elaborate the renormalization process of entropy of a nonextremal and 
an extremal Reissner-Nordstrom black hole by using the Pauli-Villars regu- 
larization method, in which the regulator fields obey either the Bose-Einstein 
or Fermi-Dirac distribution depending on their spin-statistics. The black hole 
entropy involves only two renormalization constants. We also discuss the 
entropy and temperature of the extremal black hole. 
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I. INTRODUCTION 



Bekenstein [jT| and Hawking proposed an intrinsic entropy and temperature of black 
holes. The black hole entropy and temperature are purely geometric in that the former is 
proportional to the area of event horizon of each black hole and the latter is determined by 
the surface gravity at the event horizon. Entropy is, however, a thermodynamic quantity 
that is usually defined in terms of the classical or quantum mechanical number of states 
available to a system under consideration. The geometric nature of black hole entropy has 
thus raised the fundamental questions on the origin of black hole entropy over the last twenty 
years. 

There have been various attempts to explain the black hole entropy thermodynamically. 
't Hooft calculated the entropy of a scalar field in a black hole background and found that 
with a particular brick wall just outside the event horizon the thermodynamic entropy of 
the scalar field gives the correct geometric black hole entropy 0. Quite recently Demers, 
Lafrance, and Myers (DLM) introduced a Pauli-Villars method to regularize the ultravio- 
let divergences of a Reissner-Nordstrom (RN) black hole entropy 0. In the Pauli-Villars 
covariant regularization method, one introduces bosonic and fermionic regulator fields to 
regularize the divergences. The strong point of the method of Ref. |^ is that it does not 
necessitate a brick wall to avoid possible ultraviolet divergences and yields a one-loop renor- 
malized entropy. They obtained the thermodynamic entropy of both the nonextremal and 
extremal black holes which involves only two renormalization constants. It was further 
shown that these two constants also appeared in the one-loop effective action of the scalar 
field coupled to gravity, and can be absorbed into the renormalization of the gravitational 
constant G and the coupling constants of higher order curvatures in the effective action. 
The renormalization of black hole entropy in terms of the gravitational constant and the 
coupling constants was first discussed in ^ and also discussed in |]6HTT . 



As usual, DLM used somewhat ad hoc arbitrary spin-statistics, the Bose-Einstein distri- 
bution, even for the fermionic regulator fields. In this paper we show that with the right 
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spin-statistics for the regulator fields the Pauli-Villars covariant regularization method can 
still work to yield a one-loop renormalized black hole entropy. We explicitly calculate the 
thermodynamic entropy of both the nonextremal and extremal RN black holes. The ther- 
modynamic entropy computed with this method is shown to be of the same form as those 
in [Q, but the two renormalization constants involved differ from those in |Q. 

Throughout this paper we adopt the units, c=G = h = k=l. The spacetime signature 
is (-,+,+,+). 

II. THE ENTROPY OF REISSNER-NORDSTROM BLACK HOLE 

The thermodynamic entropy of a massive scalar field in an RN black hole background 
will be studied in this section. We shall follow mostly the method used by 't Hooft . The 
RN black hole has the metric 

ds' = -(l--) (l-!±)dt' + -, ^ ^ + rW, (1) 



where r± = M ± v^M^ — Q^, and r_|_ > r_ for the nonextremal RN black hole and = r_ 
for the extremal RN black hole. The limiting case of r_ = and r+ = 2M corresponds to 
the Schwarzschild black hole. 

We consider a massive scalar field minimally coupled to the RN black hole background. 
The scalar field satisfies the (quantum) Klein-Gordon equation 
1 d'^ 1 d . d 



1/1 d , . , 



0(x) = 0. (2) 



r^^sin6'56' 89 dcjP' - 
Since the RN black hole has a timelike Killing vector field and the spherical symmetry, we 
expand the wave function by spherical harmonics, ^(x) = e~^^^Yirn{6 , (j))f{r). The remaining 
radial wave function of the Klein-Gordon equation satisfies the one-dimensional Schrodinger- 
like equation 



l d ^d. r^E^ + ^ 
— — f(r-r_)(r-r+) — j +- m 



/(r) = 0. (3) 
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The radial momentum, Pr = jr, can be found from the WKB wave function, /(r) = e''^*^^^ 



Pr 



m 



(4) 



The scalar field suffers an infinite gravitational redshift when it passes from the event horizon 
to infinity. 

We find the Helmholtz free energy from the number of states available and the statistical 
distribution. For a given E we count the action in the unit of Planck constant^ in order to 
get the number of states summed over the angular momentum states 



9iE) 



— / dmi + 1) / dr-. ^, 



X 



E'- I- — 



/(/ + !) 



1-^ 
r 



+ m 



1/2 



(5) 



A. Ultraviolet Divergence 



The radial momentum goes to infinity at the event horizon r = and it causes an 
ultraviolet divergence. To avoid the ultraviolet divergences 't Hooft introduced a brick wall 
of thickness h just outside the event horizon to regularize the divergences. The number of 
states after performing the angular momentum integration amounts to 



9{E) 



dr- 



3tT Jr++h ^1 _ ^1 _ I+y 

and the Helmholtz free energy is 



E' 



1 



1 3/2 



m 



(6) 



F = lEln(l-e-^^-). 



(7) 



For a continuous E, we may take a continuum limit to compute the Helmholtz free energy 



^Here we adopted the convention of action / drpr and the number of states 2^ / drpr as in Ref. 
PI- 
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dEg{E)- 



1 



e/3i? _ 1 • 

Substituting into (|^) and changing the variable to x = 1 — we finally obtain 



(8) 



RN 



dE 



dx- 



1 



3/2 



E — m x{l — u + ux) (9) 



and u = 1 corresponds to the exremal RN black hole. 



where m = — , e — , , , 

We focus on the possible ultraviolet divergences coming from the event horizon and 
expand the rational function part of integral (|^) around x = 0. The Helmholtz free energy 
for the nonextremal black hole is then given by 

.3 oo 



n.ext 



3n 



n=0 



dE- 



dxx 



-2+n 



E ~ m x(l — u + ux] 



3/2 



(10) 



where 



Cn 



1 ^ 2(2 - 3n) 



:i-u 



'l-u) 



(11) 



and for the extremal black hole by 



„3 oo 



ext 



dE 



3tt 



n=0 



dx X 



-4+n 



j-iz 2 2 

E — m x 



3/2 



(12) 



where 



4 ■ 5 ■ ■ ■ (3 + n) 



(n > 1). 



(13) 



We show explicitly how the ultraviolet divergent terms in the nonextremal black hole 
emerge by performing integration of the first two terms in ( PUD . By keeping only the divergent 
terms |]T^, we have 



n.ext 



37T 



71 



1 



+ 



■Infe^ 



71^(2 - 3u) 



15(1 - ii)2/34 e 8{l-u)(3^ ' ' 15(1 -u)^P^ 
where we used the integral formula, for n = 0, 1 of (p!0D , 

1 



ln(e^ 



(14) 



dx- 



X 



1 . j JE + a+x — JE + a-X 
In 



E'^ -'m^x{\-u^ux) E \^E + a+x + ^E + a-x) ' 

1 — u)m^ ± \J{1 — uYm'^ + Aum?E'^ 
2E ■ 



(15) 



We note that Fn.ext contains only a linearly divergent term and logarithmically divergent 
terms In(e^) and Ern^ In(e^) as e — 0. The Helmholtz free energy of the scalar field thus 
diverges as the brick wall thickness h (here e) collapses to zero at the event horizon. 

The ultraviolet divergent terms in the extremal black hole similarly come from the first 
four terms in ([l^ . Keeping only the divergent terms we get 



TT 



+ 



1 27rn vrW 1 2Tr^ 



+ 



ln(e2) + 



2/32 



ln(e^ 



(16) 



where the use has been made of the following integral formula [13] , for n = 1, 3 of (|T2| 



dx- 



1 



'E-VE^-m^ 



771,2x2 



X^/W 

The divergent terms are — 



2E \E + VE-^ -m^x\ 



^3ln(e2), and EmHn{e^). 



(17) 



The main purpose of this paper is to see whether one may remove these divergent terms 
and may renormalize the free energy by introducing bosonic and fermionic regulator fields 
that obey the proper Bose-Einstein or Fermi-Dirac statistics. This point will be discussed 
in Sec. III. 



B. Brick Wall Model 

't Hooft introduced a brick wall just outside the event horizon to cut off the ultraviolet 
divergence and suggested that the thermodynamic entropy of scalar fields might give the 
correct geometric entropy of black hole if a particular brick wall thickness is chosen. We 
reconsider his argument using the entropy of the nonextremal RN black hole including the 
logarithmically divergent terms. 

First, the Schwarzchild black hole case is recovered from the limiting case of r_ = and 
r+ = 2M. Repeating the calculation, focusing only on the dominant divergent terms of the 
Helmholtz free energy, and using S" = /3^|^, we find the thermodynamic entropy 



^3 

Ssch ~ — 

3% 



813 ■ 

(18) 



47r4 1 7r2m2 , , , IOtt^ , ^ 
+ — ^ln(e) - — — ln(e 



15/33 e 2(3 15/3' 
't Hooft ascribed the black hole entropy only to the linearly divergent term and prescribed 



where A = IGvrM^ is the area of the event horizon of Schwarzchild black hole. He thus 
obtained ~ j. However, assigning (3 = SttM one gets the on-shell black hole entropy 

11 m^y4 , , , 1 , , , , . 

Ssch + n e n e . 20) 

^^'^ 720 e 967r ^ ^ 180 ^ ^ ^ ^ 



It is to be noted that the second term of (|T8|) gives rise to an entropy proportional to the 
area of the event horizon, so this term should not be neglected even for small e. We should 
instead prescribe 



efl-^ln(e)Ul^ii^ (21) 
V 24% ^ ' 180A ^ ' 



to get the black hole entropy 



Ssch ~ (22) 



Similarly, we find the thermodynamic entropy of the nonextremal RN black hole 



Sn.ext 7r~ 



4%^ 1 vr^m^ , , , 87t\2 - 3m) , , / 
+ TTT. 7T75 ln(e) - ,,,, ^ ln(e) 



15(1 - u^p^ e 2{l-u) 13 15(1-^)3/33 



(23) 



The on-shell entropy, /3 = ^^^^ is 



1 — -u 1 rn^A^ , , 2 — S-u , , , 

Sn ext + n(e) n(e) . (24) 

n.ext e 967r ^ ^ 360 ^ ^ ^ ^ 



By prescribing the brick wall condition 



V 2A-K ^ ' 180v4 ' ^ ' 



we obtain the black hole entropy 



Sn.ext ~ ^- (26) 



Thus we have seen that the same entropy formula is obtained by including the logarith- 
mically divergent term that 't Hooft neglected. The contribution of logarithmically divergent 
term to the 2D black hole entropy was discussed in 



III. PAULI-VILLARS REGULARIZATION 



It is our observation that by using the correct Bose-Einstein and Fermi-Dirac statistics 
for the regulator fields we can remove the same divergent terms for both the nonextremal 
and extremal RN black hole, because they contribute to the Helmholtz free energy with 
opposite signs, and the Bose-Einstein distribution 

and the Fermi-Dirac distribution 

Jo "^^^^1 = ^^-^ 
differ by only a constant factor. This suggests that in order to remove the ultraviolet 
divergences we may introduce bosonic and fermionic regulator fields that satisfy the Klein- 
Gordon equation with arbitrarily large masses mj. The rrii will be determined later from 
the mass conditions that make all the ultraviolet divergent terms cancel each other. 

Repeating the same calculation as for the bosonic massive scalar field, we get the number 
of states for each regulator field 

g,{E) = —( dr _M ^^"^ 



(29) 



Svr Jr++h - — ) (l - — ^ 

where i labels the species of regulator fields. Similarly, we obtain the Helmholtz free energy 
contribution of the zth regulator field 

^^ = 4Eln(lTe-'^^-), (30) 

and for continuous E 

■^3 rco 1 /■ 1 r l3/2 

Fi = =1=— / dE—T^ / dx^rz — E'^ — x(l — u + ux)m^ , (31) 

where the upper (-) sign is used for bosonic fields and the lower (-I-) sign for fermionic fields, 
and we used the same variables and parameters as in Sec. II. As before, u = 1 corresponds 
to the extremal black hole case. 
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A. Nonextremal RN Black Hole 



The ultraviolet divergent terms of each regulator fields in the nonextremal RN black hole 
come from the singular behavior of the Helmholtz free energy (^TJ) near x = 0. We compute 
the Helmholtz free energy around x = 0: 

T"^ ^ \ f r 1 3/2 

Fi = T— y Cn / dE— / - mjxil -u + ux) (32) 

Btt Jo e^^T^Je L '\ 

where 



Co — y, r^, Ci — -J- (33) 

'1 — uY 1 — uY 



Regardless of spin-statistics of regulator fields, the free energy contains a linearly divergent 
term ^ and logarithmically divergent terms i?^ln(e^) and Em\\n.{e^) as e — as shown 
in Sec. II. We may remove these divergent terms by introducing the correct bosonic and 



fermionic regulator fields that contribute to the free energy (^) with opposite signs, and 
thereby regularize the free energy and entropy. 

The linearly divergent contribution to the free energy is 

where Nb and Np are the number of the bosonic and fermionic fields, respectively. We 
can remove this term by introducing 7k {k a positive integer) bosonic fields (including the 
original scalar field) and 8A; fermionic regulator fields. In this paper we shall take the 
minimum number of regulator fields, that is, 7 bosonic and 8 fermionic fields. 
Next, we consider the logarithmically divergent contributions 

r 7 \ 7r3(2-3^)r^ 



If we further impose a condition on the masses of the regulator fields 



E^?-^E"^HO, (36) 



then both the hnearly divergent term and logarithmically divergent terms are made vanish. 
So we may remove the brick wall, e = 0, under the condition (|36D, and instead regularize 
the free energy by the large masses, mj, of regulator fields. 

We find the Helmholtz free energy in the limit e = 0. For this purpose we expand the 
right hand side of the integral (]TB|), use the mass condition ( PBD to remove In(e^), and obtain 
the only nonvanishing integral 



Then the Helmholtz free energy is given by 

7 ^, 



^ n.ext T, 

611 



7 .nr, Elni- '- 

c 7-11 / (l-u)'^m'^+4um'^E'^ 
3 A „ /-oo E\n{^ '—rjrrr^ 



4(1 -u)^ ' Jo 
„ Jo 



B 

2 /'^^ZT^^^^^V 16^ 



4(1 -u)y ' Jo e^^ + 1 



:i - uY V ^0 e/^"" + 1 



(38) 



For large masses of regulator fields, we expand the logarithmic terms in (^) in the inverse 
power of masses. The leading terms are (ln(m^ — 21n(ii^)). Keeping only the surviving 
divergent contributions to the Helmholtz free energy as the regulator masses go to infinity 
we obtain 

TirlB 2-K^r%{2-?,u)A 



n.ext 



12(l-n)/32 45(l-n)3/34 ' 



(39) 



where 



^=-Eln(m?) + ^ElnK^), (40) 

B ^ F 

7 ^ 8 

B F 



TT^ ^ Jo - I Tx^ p Jo e* + 1 
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From the inequality t > ln(t) > ^ for 1 > a > 0, we may find tlie bound for the integrals 



of (|TD: 



r(3)c(3)>r*^M!^> 

JO e' — 1 



tln(t) ^ r(2-a)C(2-a) 



ea 



5r,3,c(3,>f.^.-(.-^)I(i.ifi^, 

The black hole entropy is then 



_ ^r;g SttV^ (2 - 3^)^ 



Note that if we had used the same spin-statistics even for the bosonic and fermionic regulator 
fields, the constants A and B would become 

3 3 3 3 

A= -Y.Hmf) + T.Hml), B = J2m^iln{m^i)-J2m^ilnim^i) (44) 

B F B F 

for three bosonic and fermionic fields, repectively, which are the same as those in |Q. 

The entropy obtained so far is off-shell and we substitute the Hawking temperature, 
(3 = Yz^^ nonextremal RN black hole to get the on-shell entropy 

_ B A (2 - 3u)A 
bn.ext - 4 + 180 ' 

where A = Anr'i^ is the surface area of the event horizon. It should be remarked that 
we obtained the exactly same form of entropy as in [Q, but with different renormalization 
constants A and B. A and B might be related with the renormalization of G and the 
coefficients of the effective action 

B. Extremal RN Black Hole 

Now we turn to the exremal RN black hole, the case with u = 1. The Helmholtz free 



energy contribution of the ith regulator field is 



rp-i POO 1 



/ dE— — - / dx—- r- - mfx^ . (46) 



3/2 
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We expand (^iJ^^-^i around x = 0, and perform the integration 



^3 ^ ^ ^oo 1 /. _ ,,„r_o o nl3/2 



The divergent terms are i?^ln(e^), and Emf\n{e^). Taking the correct spin- 

statistics of the bosonic and fermionic regulator fields into account, the ultraviolet diver- 
gences can be removed with the same condition on the number of regulator fields and the 
mass condition ( PBf ) as in the nonextremal RN black hole. We remove again the brick wall 
outside the event horizon, and let e = 0. Then the nonvanishing contribution to the free 
energy from the lower limit of integration is 

^e.t 9^4 ' 

and the entropy is found to be 

^^^^ 

The black hole entropy (^9|) has also the same form as in ||^, but differs only by the renor- 
malization constants A and B. 



IV. CONCLUSION AND DISCUSSION 

In this paper, we obtained the thermodynamic entropy of the RN black holes using the 
Pauli-Villars regularization method. The primary difference of our method from that of 
DLM is that we used the Bose- Einstein statistics for the bosonic regulator fields and the 
Fermi-Dirac statistics for the fermionic regulator fields, whereas they used the same Bose- 
Einstein statistics for both the bosonic and fermionic regulator fields. We confirm that the 
thermodynamic entropy of either the nonextremal or extremal RN black hole involves only 
two constants A and i3 as in [Q. But the regulator masses satisfy the mass condition ( |36D 
which is different from that in |^ due to the different spin-statistics of regulator fields. The 
renormalization constants of the thermodynamic entropy might be related with the renor- 
malization of the gravitational constant G and coupling constants of higher order curvatures 
of the effective action. 
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The thermodynamic entropy of the extremal black hole has the same form as that of the 
nonextremal black hole. There still remains the problem of defining the temperature and 
the entropy. From our result we can infer several possibilities. First, as argued in many 
literatures the entropy in (pQ]) goes to zero as the temperature approaches to zero. This fact 
is consistent with the argument in many literatures that both temperature and entropy of 
the extremal black hole must vanish. The second possibility is that the temperature is still 
inversely proportional to r+ as for the Schwarzschild or nonextremal RN black holes and the 
area rule of black hole entropy remains valid. The third possibility is that the temperature 
is arbitrary [|14| but it is related to the entropy as (^). 

In m it was shown that the renormalization constants appearing in the entropy coincide 
with those from the coupling constants of the higher order curvatures in the one-loop effective 
action. It would be interesting to find the way to relate these renormalization constants with 
the correct spin-statistics for the regulator fields. 
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